Exact 2+1 Flavour RHMC Simulations

Traditionally, “2+1” flavour simulations have been performed using the R algorithm [1], but an exact algorithm is clearly desirable. There also exist exact Polynomial Hybrid Monte Carlo algorithms [2,3], which can be used for 2+1 simulations: however, such algorithms are expensive with regard to memory consumption and for the case of asqtad fermions such an algorithm would be impractical due to the very expensive force term calculation. We explore the use of the Rational Hybrid Monte Carlo (RHMC) algorithm applied to asqtad and domain wall 2+1 simulations.